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According to the favour models for the formation of large-scale structure in 
the Universe (in which the dynamics of the Universe is dominated by cold 
dark matter), the distribution of galaxies and clusters of galaxies should be 
random on large scales. It therefore came as a surprise when a periodicity 
was reported^ in the distribution of high-density regions of galaxies in the 
direction of Galactic poles, although the appearent lack of periodicity in other 
directions led to the initial report being regarded as a statistical anomaly^. A 
subsequent study^“® also claimed evidence for periodicity on the same scale, 
but the statistical significance of this result was uncertain due to small number 
of clusters used. Here, using a new compilation^ of available data on galaxy 
clusters, we present evidence for a quasiregular three-dimensional network of 
rich superclusters and voids, with the regions of high density separated by ~ 
120 Mpc. If this refiects the distribution of all matter (luminous and dark), 
then there must exists some hithero unknown process that produces regular 
structure on large scales. 

During the past few years the number of clusters with measured redshifts has increased 
considerably. To search for the possible presence of a regularity of the distribution of matter 
in the Universe we have used a new compilation^ of available data on rich clusters of galaxies 
catalogued by Abell and collaborators®’®. The compilation has made use of all (~300) 
published references on redshifts of both individual galaxies and Abell galaxy clusters. 
Individual galaxies were associated with a given Abell cluster if they lay within a projected 
distance of <1.5 h~^ Mpc (1 Abell radius) and within a factor of two of the redshift 
estimated from the brightness of the cluster’s 10-th brightest galaxy, using the photometric 
estimate of Peacock & West^® {h is the Hubble constant in units of 100 km/s/Mpc). The 
compilation contains measured redshifts for 869 of the 1304 clusters with an estimated 
redshift up to 2 ; = 0.12. For the present analysis we used ah rich clusters (richness class 
i? > 0) in this compilation with at least two galaxy redshifts measured. The omission of the 
435 clusters without measured redshifts does not affect our result because an appropriate 
selection function was used. 

This cluster sample (including clusters with estimated redshifts) was used to con¬ 
struct a new catalogue of 220 superclusters of galaxies^^’^^. These are systems of clusters 
where the distances between nearest neighbours among member clusters do not exceed 
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FIG. 1. The distribution of 319 clusters in 25 very rich superclusters with at least 8 mem¬ 
bers (including 58 clusters with photometric distance estimates), illustrating the network 
in the cluster distribution in supergalactic coordinates. In the lower panels clusters in 
the northern and southern Galactic hemispheres are plotted separately. The supergalactic 
y = 0 plane coincides almost exactly with the Galactic equatorial plane, that is, with the 
zone of avoidance due to galactic absorption. The grid with step size 120 h~^ Mpc corre¬ 
sponds approximately to distances between high-density regions across voids. In the two 
upper panels and in the lower right panel several superclusters overlap due to projection 
but are actually well-separated in space. 

24 h~^ Mpc. In this way we hnd high-density regions in the distribution of clusters of 
galaxies. In Fig. 1 we plot clusters located in rich superclusters with at least 8 member 
clusters. We see a moderately regular network of superclusters and voids with a step size of 
^ 120 ±20 h~^ Mpc where chains of superclusters are separated by voids of almost equal 
size^^d2 whole distribution resembles a three-dimensional chessboard^^. Nearest 
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neighbour test, and pencil-beam and void analysis indicate that clusters in poor super¬ 
clusters with less than 8 members and isolated clusters form a more uniform population, 
preferentially located in void walls between rich superclusters but not hlling the voids^^. 

To quantify the regularity of the cluster distribution we have calculated the correlation 
function and the power spectrum for clusters of galaxies. The correlation function describes 
the distribution of clusters in the real space, the power spectrum in the Fourier space of 
density waves. Analysis of various geometric models has shown that if superclusters form a 
quasiregular lattice with an almost constant step size then the cluster correlation function 
is oscillating, it has alternate secondary maxima and minima, separated by half the period 
of oscillations. The period of spatial oscillations of the correlation function is equal to the 
step size of the distribution^^. The amplitude of the power spectrum at the wavelength 
corresponding to that period is enhanced with respect to other wavelengths, that is, it is 
peaked. In contrast, if superclusters are located randomly in space then the correlation 
function approaches zero level at large separations and the power spectrum turns smoothly 
from the region with positive spectral index on large wavelengths to a negative index on 
small wavelengths^^. 

To calculate the power spectrum we have used the sample which contains clusters 
with measured redshifts only and lying in both Galactic hemispheres out to the distance 
covered by our cluster and supercluster catalogues. The power spectrum was derived 
using two different methods, a direct one where we calculate the distribution of clusters 
in the wavenumber space, and an indirect method where we hrst calculate the correlation 
function of clusters of galaxies and then hud the spectrum. In the latter case we make use 
of the fact that the power spectrum and the correlation function are related by the Fourier 
transform. 

In both methods the main problem is the calculation of the selection function of 
clusters of galaxies which corrects for incompleteness both at low Galactic latitude b and 
at large distances r from the observer. The selection f un ction can be represented by linear 
functions of sin 6 and r (for numerical values of selection f un ction parameters see refs 12 
and 15). Both methods to derive the power spectrum yield similar results. However, 
parameters of the correlation function are insensitive to small inaccuracies of the selection 
function^®, thus the indirect method yields more accurate results for the power spectrum. 

The power spectrum for clusters of galaxies is shown in Fig. 2. This power spectrum 
is among the hrst for three-dimensional data having measurements on scales well above 
100 h~^ Mpc. On very large scales the errors are large due to incomplete data. On moder¬ 
ate scales we see one single well dehned peak at a wavenumber ko = 0.052 h Mpc~^. Errors 
are small near the peak, and the relative amplitude and position of the peak are determined 
quite accurately. The wavelength of the peak is Xq = 2%/ko = 120 ±15 h~^ Mpc. Near the 
peak, there is an excess in the amplitude of the observed power spectrum over that of the 
cold dark matter (GDM) model (see below) by a factor of 1.4. Within observational errors 
our power spectrum on large scales is compatible with the Harrison-Zeldovich spectrum 
with constant power index n = 1, and on small scales with a spectrum of constant negative 
power index n = —1.8. 

The power spectrum is often used to compare the distribution of matter in the Universe 
with theoretical predictions. Gurrently popular structure formation theories based on the 


3 




0.01 0.10 

k [h Mpc^'] 

FIG. 2. The power spectrum 
P{k) for 869 clusters with measured redshifts is plotted with solid circles. The spectrum is 
calculated from the cluster correlation function via the Fourier transform, both functions 
are found in redshift space. Errors were determined from 2a errors of the correlation func¬ 
tion found from the scatter of different simulations. As for each wavenumber we integrate 
over the whole interval of the correlation function individual values of the spectrum are 
not independent of each other. To check the indirect method of calculation of the power 
spectrum we have used simulated cluster samples having similar selection effects as real 
ones. This check was performed for a wide variety of models with different initial spectra. 
Our results show that the true spectrum can be restored over the wavenumber interval 
from k ~ 0.03 h Mpc“^ towards shorter waves until k ~ 0.3 h Mpc“^. The solid line is 
the standard CDM {h = 0.5, 0 = 1) power spectrum enhanced by a bias factor of 6 = 3 
over the four year COBE normalisation. 

dynamics of a Universe dominated by CDM yield spectra of which one example is plotted 
in Figure 2 as a solid line. The spectrum is rising on long wavelengths A (small values 
of the wavenumber k = 27r/A), and falling at short wavelengths (large values of k). The 
transition between short- and long-wavelength regions in the CDM-spectrum is smooth. 
The distribution of superclusters in CDM-models is irregular^®. 

As we see, the relative amplitude of the observed power spectrum above the standard 
CDM-type model is not very large. Thus we may ask the question: within the framework 
of the standard cosmogony, how frequently can we expect to hnd a distribution of clus¬ 
ters which has a power spectrum similar to that observed? To answer this question we 
determined the correlation function and power spectrum for clusters in rich superclusters 
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of CDM-type models. In the spectral range of interest the power spectrnm of the standard 
CDM model is similar to the spectrnm of a random snperclnster modeh^. We make nse of 
this similarity by generating 1,000 realizations of the random snperclnster model, applying 
the selection effects as fonnd in clnster distribntion, and determining the parameters of the 
clnster correlation f un ction and power spectrnm. We measnre the mean period and ampli- 
tnde of the oscillating correlation fnnction and their respective scatter. We also calcnlate 
the deviations for individnal periods. This test shows that combination of parameters (for 
instance period vs. amplitnde) close to the observed valnes occnrs in approximately 1 % 
of cases, bnt the simnltaneons concnrrence of all parameters with observations is a mnch 
more rare event. Thns some change in the initial spectrnm of matter is necessary in order 
to explain the observed correlation fnnction and power spectrnm for clnsters of galaxies. 

The regnlarity of the distribntion of snperclnsters is qnite striking and has some 
similarity with the regnlarity fonnd in pencil-beam galaxy snrveys^. Independent evidence 
for the presence of a preferred scale in the Universe aronnd 100 h~^ Mpc comes from a 
core-sampling analysis of the distribntion of sheet-like and hlamentary strnctnres^^, and 
from a power spectrnm analysis of the two-dimensional distribntion of galaxies^® in the Las 
Campanas Redshift Snrvey. Available data, however, are insnfhcient to say whether one-, 
two-, and onr three-dimensional snrveys measnre identical physical scales in the Universe, 
or, if so, whether differences in n um erical valnes of the scale are dne to differences in 
methods applied in determination or to a real cosmic variance in the scale. 

According to onr present nnderstanding of the Universe the observed distribntion of 
Inminons matter is strongly correlated with the density pertnrbations of some dark non- 
baryonic matter component(s) at the moment of recombination. These pertnrbations are 
dne to processes in the very early Universe so that their power spectrnm at recombination 
depends both on these processes and the fnrther evolntion of the Universe, in particnlar 
the transition from radiation to matter dominated expansion. We have shown that the In¬ 
minons matter in galaxy clnsters is mnch more regnlarly distribnted than expected. Thns 
we end np with the conclnsion that onr present nnderstanding of strnctnre formation on 
very large scales needs revision. 
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